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^ ■ ABSTRACT 

. This paper aims to connect the theory of relativistic cosmology number counts 

with the astronomical data, practice, and theory behind the galaxy luminosity 
function (LF). We treat galaxies as the building blocks of the Universe, but 
^ . ignore most aspects of their internal structures by considering them as point 

0\ . sources. However, we do consider general morphological types in order to use 

^ . data from galaxy redshift surveys, where some kind of morphological classifica- 

O , tion is adopted. We start with a general relativistic treatment for a general space- 

Q . time, not just for Friedmann-Lemaitre- Robertson- Walker, of number counts, and 

then link the derived expressions with the LF definition adopted in observational 
cosmology. Then equations for differential number counts, the related relativis- 



2 ■ tic density per source, and observed and total relativistic energy densities of the 



c/3 , universe, and other related quantities are written in terms of the luminosity and 

^ ' selection functions. As an example of how these theoretical/observational rela- 



X 
J3 



tionships can be used, we apply them to test the LF parameters determined from 
the CN0C2 galaxy redshift survey, for consistency with the Einstein-de Sitter 
(EdS) cosmology, which they assume, for intermediate redshifts. We conclude 
that there is a general consistency for the tests we carried out, namely both 
the observed relativistic mass-energy density, and the observed relativistic mass- 
energy density per source, which is equivalent to differential number counts, in 
an EdS Universe. In addition, we find clear evidence of a large amount of hidden 
mass, as has been obvious from many earlier investigations. At the same time, 
we find that the CN0C2 LF give differential galaxy counts somewhat above the 
EdS predictions, indicating that this survey observes more galaxies at 0.1^2:^0.4 
than the model's predictions. 

Subject headings: cosmology: number counts, luminosity function, relativity 
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1. Introduction 

In cosmology, it is one thing to obtain universe models from solutions of Einstein's 
field equations, and quite another to express these solutions in terms of quantities and 
parameters which can be related directly to real astronomical observations. This intrinsic 
difficulty of cosmological research was perceived long ago by its practitioners, more or less 
by the time modern cosmology was launching its own foundations, in the 1920's and early 
1930's. This perception initiated a process which led to the estabhshment of a division 
between theoreticians and observers, with the former being grouped in what may be termed 
as relativistic cosmology, usually considered to be the realm of the first task above, and the 
latter forming observational cosmology, which lies in the realm of the second one. Despite 
this division, there has always been a few, of course, who were able to link successfully 
these two sub-fields {e.g., Tolman 1934; McCrea 1934, 1939; Sandage 1961, 1995; Elhs 1971; 
Weinberg 1972; Elhs and Perry 1979; Ellis et al. 1984; Peebles 1980, 1993; Longair 1995). 
And more recently, there have been many individuals and groups of researchers who have 
successfully and closely hnked cosmological observations (particularly those of the cosmic 
microwave background radiation and those involving deep galaxy surveys) to the standard 
Friedmann-Lemaitre-Robertson- Walker (FLRW) theory. However, more needs to be done, 
particularly with regard to the broader issues of using observations to determine how good 
a model FLRW is on various cosmological length scales. 

Each of these two approaches to cosmology (relativistic-theoretical and observational) 
has its own distinctive features. Relativistic cosmologists have been so successful in coping 
with the difficulties of obtaining exact solutions of Einstein's field equations, that today 
one of the problems of this area is not the lack, but rather the excess, of solutions without 
the correspondent astrophysical modeling (Krasihski 1997). So, in relativistic cosmology 
a different cosmological model means in fact a different spacetime metric. Observational 
cosmologists, on the other hand, usually work almost exclusively with FLRW spacetimes, 
since, by its very nature, observational verification in cosmology is a formidable task, and 
since there are indications that the Universe is nearly FLRW on the largest scales. The fact 
that the characterization of FLRW geometry can be reduced to a set of a few parameters to be 
determined observationally, has allowed observational verification in cosmology to be reduced 
to manageable proportions. Or so it has been thought. Consequently, for observational 
cosmology a different cosmology often means, in fact, different parameter values for an 
FLRW spacetime metric, changing, for instance, its general topology from fiat, to closed or 
open. Nevertheless, however different, these two sub-areas share the same general purpose of 
modern cosmology, which is to determine the spacetime geometry and the mass distribution 
of the Universe. 



-3- 



One of the most recent and comprehensive attempts to unify these two approaches in 
cosmology has been, perhaps, the ideal observational cosmology program, which proposes 
to characterize in detail the way in which cosmological observations can be directly used 
to determine the cosmological spacetime geometry (Ellis et al. 1985; Nel 1987). The main 
motivation behind this program has been to determine the spacetime metric of our universe 
directly from astronomical observations, without assuming a cosmological model beforehand 
(see EUis, Matravers and Stoeger 1995). In the process of doing that, one must first determine 
what is and what is not decidable in cosmology on the basis of astronomical observations 
(Ellis et al. 1985, p. 317). 

However, despite its obvious appeal, the ideal observational cosmology program still 
encounters a very basic problem: the astrophysical evolution, both in number and in lu- 
minosity, of galaxies, which are the fundamental astronomical makers of cosmic spacetime. 
Astrophysical evolution and cosmology are convolved in every presently practical astronom- 
ical observation of cosmological relevance. In order to implement the ideal observational 
cosmology program, not only do we require precise observations of redshift, galaxy number 
counts, and observer area distances (or luminosity distances), in the spherically symmetric 
case, but we also need an adequate model of how the populations and luminosities of the 
different types of galaxies evolve with redshift. It is presently impossible to determine this 
without assuming some cosmological model. This problem is so difficult that astronomers 
usually assume that the universe is well described by the standard FLRW geometry, and 
analyze, interpret and present their observational results in terms of that class of models, 
with or without a cosmological constant. On the largest scales - that is, on the largest scales 
probed by the cosmic microwave background radiation (CMBR) - there is good reason for 
doing this (see Stoeger, Maartens and Ellis 1995), although on scales less than these, this 
assumption is still in need of more adequate justification, inasmuch as we are incapable so far 
of determining the smallest length scale on which the FLRW approximation is valid. Thus, 
most observational results in cosmology either explicitly or implicitly assume an FLRW 
model in their presentation. 

Although this approach has enabled the rapid progress of cosmology to this point, it 
may significantly hamper taking the next step: determining observationally how significant 
the deviations from FLRW are on various length scales, once we have some adequate model 
of galaxy evolution. 

Finally, in some cases, in order to obtain simplified models, the approach taken in 
observational cosmology - particularly at low redshifts - has not been fully relativistic, and 
sometimes involves Euclidean simplifications, without proper verification of the limits within 
which such simplifications are justifiable (Ellis 1987; Ribeiro 2001). Though with the advent 
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of much deeper surveys this is very rare now, it still causes some confusion, since some 
standard references show simplified presentations which, often for clarity, do not include 
these corrections.^ 

The aim of this paper is to bridge the relativistic and observational cosmology ap- 
proaches. But, it is much less ambitious than the ideal observational cosmology program, 
although we do intend to start with a fully relativistic approach and go as far as possible with 
the general equations, only then specializing them to particular geometries. In this paper, as 
an important example, we specialize to the Einstcin-dc Sitter model. The basic philosophy 
follows the general proposal for limited bandwidth observations of cosmological point sources 
outlined in Ribciro (2002). However, our focus here will be on galaxy number counts, since 
this a very important cosmologically relevant quantity, giving information about the density 
of mass-energy in the universe, and can be determined by careful analysis of the data in 
galaxy redshift surveys. More particularly, therefore, our purpose is to present the detailed 
relationship between galaxy number counts, or galaxy luminosity functions, and relativistic 
mass-energy density and relativistic mass-energy density per source of the universe for any 
general cosmological model, and, as an example of how these relationships can be used, 
specialize them to FLRW and employ them to test the consistency of the FLRW model, in 
this case Einstein-de Sitter (EdS), assumed in deriving luminosity functions from a given 
redshift survey (the CN0C2 intermediate redshift survey) with the mass-energy density and 
mass-energy density per source, as functions of redshift, implied by that model. Do the 
number counts really determine these parameters to be EdS, as they should? 

There are other ways, of course, in which we can use this framework. For instance, 
investigating the interplay of evolution, selection-effect and incompleteness models with var- 
ious cosmologies in interpreting the same data, and testing simplifying and non-relativistic 
assumptions for certain redshift regimes. It should be mentioned here, too, that the other im- 
portant cosmologically relevant observational parameter which complements number counts 
in determining the cosmological and evolution models is observer area distance (or luminos- 
ity distance). An scheme analogous to the one we present here for number counts can be 
established for this parameter (see Stoeger 1987; Stoeger and Araiijo 1999). 

Due to the reasons stated above, we shall in this paper specialize our model to FLRW. 



^ A common Euclidean approximation used in the past in cosmology is to calciilate the galaxy number 
counts by means of an energy flux density which is not corrected by the redshift (see, e.g., Peebles 1993, p. 
214). Another commonly used Euclidean approximation was to take the mass density p to be constant when 
z < 0.1 without a proper verification if that is theoretically advisable within the problem under discussion. 
Although such approximations, often used in the past, are a rarity in recent redshift surveys analysis, the 
fact that they can be found in often-quoted texts is still a source of confusion. 
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But we intend to expand this study eventually to non-FLRW metrics. This is justified 
because we wish to see what are the main differences between the predictions of reasonable 
non-FLRW metrics and those of the standard (FLRW) model. The only way we can really 
determine how good the Friedmann models arc, is to fit almost-FLRW and non-FLRW 
models to the data, and compare them with the Friedmann fits to see if these latter are 
better (Ribeiro 2002). As the next logical step the obvious choice would be the Lemaitre- 
Tolman-Bondi spacetime, since this is the most widely used metric for testing cosmological 
models after the standard model (Krasinski 1997). 

Much of this detailed theoretical relativistic treatment has been presented before, but 
it was never linked or applied directly to luminosity function (LF) and number count data, 
with a careful discussion of the problems and obstacles of doing so. This is the first in a series 
of papers investigating how a fully rigorous general theoretical cosmological framework can 
effectively incorporate cosmologically relevant observational data - first assuming FLRW and 
testing for consistency, and then testing the FLRW hypothesis itself on various cosmologically 
relevant length scales by seeing if non-FLRW models also fit the data, with reasonable 
evolution, selection-effect and incompleteness assumptions. 

As mentioned above, since most cosmological data is analyzed and presented employing 
the assumption of FLRW, we shall first perform consistency tests of the derived LF's with 
the FLRW model they assume. Models with non-zero cosmological constant will be discussed 
in a follow-up paper (Stoeger and Ribeiro 2003) 

The main conclusions we reach in this study is that the luminosity function parameters 
derived from data acquired by the recent rcdshift surveys are roughly consistent with the 
theoretical predictions of the standard Einstcin-dc Sitter (EdS) model. However, the rela- 
tivistic energy density per source, or equivalently the galaxy number counts, departs from 
this model's predictions at 0.1^ 0.4. That may indicate either an inadequacy of the EdS 
cosmology to represent the data at those redshifts, or the rough number evolution model 
used in the CN0C2 LF calculations needs improvements. 

This paper is organized as follows. In §2 we derive metric independent relations based 
on number source counting and the luminosity function. In §3 we collect the available data 
on the luminosity function and organize them in a form usable to our purposes. Section 
4 specializes the general equations to the Einstein-de Sitter spacetime geometry, whereas 
section 5 carries out two consistency tests between the observational data and Einstein-de 
Sitter's cosmological model theoretical predictions. Section 6 discusses the results and their 
implications. 
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2. Model Independent Counts of Cosmological Sources 



We begin by deriving completely general cosmological number counting dependent quan- 
tities, which could, in principle, be determined directly from observations. 



The key result derived by Ellis (1971) for the number of cosmological sources in a volume 
section at a point P down the null cone is 



Here n is the number density of radiating sources per unit proper volume in a section of a 
bundle of light rays converging towards the observer and subtending a solid angle dJlo at the 
observer's position, c?^ is the area distance ^ of this section from the observer's viewpoint, 
is the observer's 4-velocity, k"" is the tangent vector along the light rays, and y is the affine 
parameter distance down the light rays constituting the bundle (see figure 1). Notice that n 
is measured in the rest-frame of the counted sources. This equation is metric independent, 
being, therefore, valid for all cosmological models (see Ellis 1971, p. 159 for more details). 
It should be noted that, in general, n has a dependence both on the cosmology and on the 
number evolution of the sources. Only when the density of galaxies is constant in a comoving 
volume will n be free of number evolution. 

While the equation above is precise, it is written in terms of unobservable quantities, 
like the affine parameter y. It can be expressed in terms of observable quantities, but this 
has been previously done in observational coordinates, which are quite different from the 
usual ones (EUis et al. 1985). 

The general definition for the redshift can be written as (EUis 1971, p. 146), 



2.1. Relativistic Density of the Universe 



2.1.1. Relativistic Number Counts 



dN = {dAfdno[n {-k''ua)]p dy. 



(1) 



1 + z 




[■""fca] observer [''^"'^a\c{q) 



(2) 



Substituting equation (2) into equation (1) yields. 




(3) 



^ In the literature is also known as angular diameter distance, observer area distance, and corrected 
luminosity distance (Ribeiro 2001, footnote 9, p. 1711). 
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AU the light rays reaching the observer's worldhne C at the same time form a hyper- 
surface which we call a past light cone, and which can be represented by a constant quantity 
over this surface. Following the notation of Ellis' et al. (1985, p. 324), let us call the function 
generating this surface w. Therefore, 

w — w{x"') — const, (4) 

defines the past hght cones along the observer's worldhne, the locus of light rays reaching 
the observer at a given moment of time. This constant is arbitrary, and we will choose it to 
be equal to the observer's proper time r at its worldhne C. So, we may write the following 
equation, 

= Mc- (5) 

If we now define a vector field ka, such that, 

then this vector field generates past light cone surfaces. The minus sign was chosen because 
we are interested in incoming light rays, which are indeed geodesies of null type if their 
generating vector field obeys the following conditions (Ellis 1971; Ellis et al. 1985; Schneider 
et al. 1992, p. 94), 

k'^K^O, A;V = 0. (7) 
Now, using the definitions of ka and u"', it follows that, 

dw dx"' dw 

Considering equation (5) we then conclude that, 

dw 



dr 



1, (9) 



c 



and 

= -1. (10) 

This result allows us to write equation (3) as below, 

dN = (d^)^(l + z)dno n{y) dy. (11) 

This equation means that the local displacement at the rest-frame of the galaxy can be 
rewritten in terms of the redshift, as follows, 

dl = [-u'^kalp dy^(l + z)dy. (12) 

We can actually write equation (11) completely in terms of 2;, with n and as functions of 
since z — z{y, 9, </?). This is an important point to bear in mind, since in the next pages 
we will often refer to the quantity dN/dz, but rarely to dN/dy. 
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2.1.2. Average Galactic Mass 

Although the number density of sources is very important for testing cosmological mod- 
els, so is the mass-energy density of the universe. Thus, for completing our comparison of 
theory with observational data, we must at some stage relate mass-energy density to the 
number density of sources. To accomplish this we need to determine the average galaxy rest 
mass. At present this is not possible to determine from observations. So, let us indicate 
how it can be determined in principle, as well as discussing briefly the rough estimate we 
use here, for the sake of illustrating our general procedure. 

Let us define A4g to be the average galaxy rest mass. Since galaxies of different morpho- 
logical types may have masses varying from 10^ to 10^^ solar masses,^ we may think of A4g as 
the simple average M.g — (i/V) X]^=i -^v, where V is the number of galaxy morphological 
types being considered (ellipticals, spirals, irregulars, etc) and M.^ is the typical rest-mass 
value for galaxies of a certain morphological class v. By "typical" we mean a simple average 
mass vahie. Note that this definition just takes the masses measured in nearby galaxies and 
extrapolates them to higher z. 

Nonetheless, galactic class population abundances may vary with the redshift due to 
galaxy evolution, that is, galaxies may change from one morphological type to another in 
different epochs. In addition, galaxy rest masses are measured in the source's rest frame. 
Furthermore, we may consider that some morphological types will be present only in certain 
redshift ranges, and be absent in others. Therefore, in general A4v will be a function of z, 
since the typical mass of a galaxy of a given type will vary with redshift. M.v{z) will still be 
an average of sorts - the average mass of, for instance, spiral galaxies in a redshift bin Az. 
What M.v{z) is would have to be determined observationally and/or astrophysically, that is, 
with a combination of both observation and theory. 

Considering this, it is therefore more appropriate to define Adghy means of the following 
weighted average, 

V 

Mg{z) = J2Pv{z)M,{z), (13) 

v=l 

where is the galactic morphology population fraction, giving the abundance of each galactic 
type relative to the total number of counted galaxies in each redshift range. This quantity 



^ These values are just a rough estimate, based on Elmegreen (1998) and Sparke and GaUagher's (2000) 
data, as dwarf galaxies may have smaller masses. 
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is normalized by 

V 

Y.Pv{zf) = l, (14) 

v=l 

where Zf stands for a fixed redsfiift value, or a certain redshift interval or bin A2;, to be 
determined observationally. 

Equation (13) as defined above does not account for or describe galaxy mergers, as the 
basic information provided by it are the proportions of each galaxy population in the sample's 
redshift range. However, it does allow for mergers when M.g{z) is considered together with 
the luminosity function (see §2.2 below). In the simplest possible case, we may assume just 
one morphological class {V = 1), that, of course, will not change with the redshift, and its 
mass will be equal to the simple average of all galactic rest masses available. In this way 
equation (13) reduces to 

Mg « 1O^^A<0. (15) 

Note again that this factor is just a rough estimate based on data presented by Sparke and 
Gallagher (2000, p. 204, 264; see their Table 5.1). 

As stated above, we have written M.y{z) as a function of the redshift since galaxy 
interactions, including mergers, will affect the average galaxy masses of a given morphological 
type as the universe evolves. Similarly, Pv{z) also depends on ^ - as galaxies may change 
from morphological type to another as we go from earlier epochs to later ones. Typically, 
we find more spiral galaxies at higher redshifts, and more ellipticals at lower redshifts, for 
instance. But how are we to determine M.v{z)l In principle, this can only be consistently 
done by using the galaxy luminosity functions for each type of galaxy as a function of redshift 
(which are defined and discussed below in §2.2), along some rehable measure of the masses 
of the galaxies in the sample - for instance the mass-to-luminosity ratio for them (see §2.4 
below). Then, one could integrate the luminosity multiplied by the mass-to-luminosity ratio 
and the galaxy luminosity function over the luminosity, which is just the luminosity density 
convolved with the mass-to-luminosity ratio (see eq. 48 below), for each morphological type 
in a given redshift interval, to obtain the mass density contributed by that sample of galaxies. 
This can then be divided by the integral over the galaxy luminosity function itself for that 
particular type in the same redshift interval, giving the number density of galaxies of that 
type, to obtain the average galaxy mass for that sample. Perhaps, we shall eventually be 
able to do this. However, at present there are several problems. First of all, the luminosity 
functions diverge a the faint end, and it is not clear how to choose reliably a faint-end cutoff. 
There are a large number of dwarf galaxies, and we do not know where they cease in the 
luminosity function. Furthermore, we cannot reliably determine mass-to-luminosity ratios 
which will apply to the full range of galaxies in the samples. 
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Thus, converting galaxy number densities to mass-energy densities presents us with a 
serious problem. We plan to investigate how this might be done with some confidence in 
later work. For this paper, in order to illustrate our general methods, we simply set the 
values of the Ai^{z) on the basis of the rough estimates for typical galaxies of given type v 
as provided by Sparke and Gallagher (2000) and Elmegreen (1998). These values are given 
in §5, where we apply our two consistency tests. This is by no means adequate, but our 
focus in this paper is to illustrate how our general framework can be used, when we have 
enough reliable observational data. 



2.1.3. Relativistic Density per Source 
The number density of all cosmological sources in their own rest frames can be written 

as, 

niy) = (16) 

where p is the local density as given by the right hand side of Einstein's field equations. 
This presumes that the cosmological constant is zero, and that there is no dark matter not 
associated with the sources we are considering. If there is, then we can simply add those 
components to the matter density of the sources to obtain p. Now, following Stoeger et al. 
(1992, §3) we shall define the relativistic density per source ix{z) in terms of the differential 
number counts dN/dz as 

( \ MgdN 
a[z) = ;— • (17) 

Thus, this quantity can, in principle, be determined observationally (Loh and Spillar 1986; 
Newman and Davis 2000). Remembering equation (11), we can also write /x(^) in terms of 
the cosmological density p, 

Mz) = (l + ^)p^. (18) 

Interestingly, pl{z), or dN/dz, can be precisely determined theoretically for any cosmological 
model {e.g., for FLRW) and then compared with what is measured observationally (see 
below). In addition, one can clearly see that it is easier to find iJ,{z) from observations by 
means of equation (17), while it is easier to find p{z) from theory by means of equation 
(18). It goes without saying, from what we have indicated above, that both n{z) and dN/dz 
depend on both cosmological evolution and astrophysical number evolution. 

The relationship between the observed differential counts and the total count may be 
written as 

(f).-(^)' 
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where J is the completeness parameter, giving the percentage of galaxies effectively observed 
in a galaxy redshift survey relative to the estimated total content of the sample. Thus, 

< J < 1, (20) 

and equation (17) becomes, 

Ma fdN\ 

Notice that this equation is bolometric, since the observed differential count on the right 
side is not restricted to any pass-band wavelength observations. We can similarly write the 
observed relativistic density of the universe ij.o{z) in terms of the observed differential count 
as follows, 

M„ fdN\ 

inasmuch as, 

fXo{z) = Jfx{z). (23) 

A final remark: although the equations above were written in terms of the area distance, 
expressions for the other observational distance definitions can be easily obtained by means 
of Etherington's (1933) reciprocity theorem, valid for any cosmological metric. It may be 
written as follows (see also EUis 1971, pp. 153-156; Schneider et al. 1992, pp. 110-114), 

d^^d^{l + zf^da{l + z), (24) 

where d^ is the usual luminosity distance, and da is the galaxy area distance.^ The area 
distance is, in principle, an observationally determinable quantity, as its definition is given 
by (Ellis 1971, p. 153), 

where dSo is the cross-sectional area at point P of a bundle of null geodesies converging 
towards the observer (see figure 1). So, if, by means of some astrophysical model, we are 
able to infer the intrinsic cross-sectional area dSo of some object, and measure its subtended 
solid angle dflo, then d^ can be found. Notice, however, that the cross-sectional area ^5*0 is 
defined in the rest-frame of the source. 



^ In the literature do is also known as angular size distance, effective distance, transverse comoving 
distance, and proper motion distance (Scott et al. 2000, p. 644; Ribeiro 2001, footnote 9, p. 1711). 
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2.2. Galctxy Luminosity Function 

The galaxy luminosity function is defined to be the /oca/ number density of a// galaxies 
with absolute luminosity between L and L + dL at redshift z (Weinberg 1972, p. 452; Peebles 
1993, p. 119; Peacock 1999, p. 399), 

'-^ = mM, (26) 

where 

£ = L/L,, (27) 

is the luminosity scale parameter, and dV is the proper volume element along the null 
cone, at some redshift, where the luminosity function is determined.^ Integration of above 
a certain luminosity threshold L defines the selection function ip (Peebles 1993, p. 214), 

POO 

^ [e{z)] = / mdi: (28) 
Je{z) 

while summing over all luminosities yields, 

poo 

^{z) = / (l){i)di. (29) 
Jo 

This simplified formulation is an idealization. There is no universal galaxy luminosity func- 
tion. In actuality, there are, instead, different luminosity functions for different types of 
galaxies and for different observational band- widths, and these evolve with redshift. Here 
we shall continue to let the unindexed luminosity function represent the distribution of all 
the galaxies, and designate the luminosity functions of certain classes of galaxies in certain 
observational band- width with subscript and superscript indices, respectively. 

The introduction of the luminosity and selection functions means that in order to be 
detected and counted a galaxy must have a minimum fiux, or apparent luminosity F, which 



® In the past there was sometimes confusion or unclarity about which volume definition was being used 
in a given survey when one studies the Hterature on galaxy luminosity function determination. The proper 
volume was initially the one mostly used, but when the surveys started to probe at redshifts deeper than 
0.1 or 0.2, most observers switched to comoving volume usage instead, which is now the prevalent volume 
definition among observers (Chris Impey 2002, 2003, private communications). Theoreticians, on the other 
hand, often adopt the proper volume definition. Therefore, in order to have a clear link between theory 
and observation, in this paper we shall adopt the proper volume as our volume definition for the himinosity 
function determination, unless stated otherwise. Wherever another volume definition becomes necessary in 
order to compare our theory with observational data, we will employ an appropriate conversion factor (see 
§4.3 below). 
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is related to its intrinsic luminosity L, at some redshift by the well-known expression 



L{z) 



F = . ..1 , (30) 



for the case of bofomctric flux. For a speciflc bandwidth range W, the apparent luminosity 
is further restricted to (Ellis 1971, p. 161; Ribeiro 2002), 



L{z) 



A7r{d^y{i + zy 



POO 

/ W{iy)J[u{l + z)]diy, (31) 
Jo 



where W{i') defines the spectral interval of the observed flux (for instance, the standard 
UBV system) and is the source spectrum function, that gives the proportion of radiation 
emitted by the source at some rest-frame frequency, but redshifted and detected by the 
observer as + z)] (Ellis 1971, p. 161; Ribeiro 2002). 

These two expressions for F take into account relativistic effects. If those effects were to 
be left out, we would have to use their Euclidean counterparts, suitable for a non-expanding 
universe approximation. In this case the bolometric flux is given by, 

F = — (32) 

which could also be found by using the reciprocity theorem (24) in equation (30). Similarly, 
the Euclidean equation for can also be found. Using the equation above for determining 
the luminosity function will strongly affect the galaxy counts, as the lower limit of the integral 
in equation (28) will be changed (Peebles 1993, p. 214). 

The incompleteness of the sample is indicated by the lower limit of the integral in 
equation (28), since faint galaxies with intrinsic luminosity below L{z) will not be counted. If 
we now suppose that the observationally determined luminosity function for brighter galaxies 
may also be used to account for those faint galaxies, we may then write the completeness 
parameter as, 

J[z) = 1 - — - / (i){i)di, (33) 







which agrees with condition (20), as it should. In addition, the selection function may be 
written as, 

^^;[i{z)]^J{z)<lf{z). (34) 



In order to provide a relativistic version for equation (26), we must link the definitions 
above to equation (11). However, it is clear from their definitions that the only difference 
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between n, as in equations (1), (3) and (11), and is that the latter restricts the rest-frame 
volume count to a luminosity range. Considering this, we can only conclude that. 



Notice that this equation assumes that both quantities above are measured relative to 
proper volume elements. If, however, ^ is being observationally determined by means of 
another volume definition, like the comoving volume, an appropriate conversion factor must 
be introduced in order for equation (35) be correct (see §4.3 below). 

In addition, the functional relationship between the redshift z and the affine parameter y 
remains unknown in the general case, and must be determined by solving the field equations 
(in particular, the null Raychaudhuri equation). Thus, once a cosmological metric is assumed 
or specified, it can be determined. In other words, to fix the function y — y{z) we require 
a cosmological solution of Einstein's field equations.^ Of course, in practice we may be 
using some null radial coordinate other than the affine parameter, say r. Then we would 
need r = r{z) or its inverse relationship. Below, in dealing with the Einstein-de Sitter 
cosmological model, we shall use an r which is closely related to the radial coordinate which 
is often used in describing FLRW models. 

We can write the relativistic number counts in terms of the luminosity function by 
substituting the result above into equation (11). The result yields. 



Hence, the observed differential number counts for sources whose intrinsic luminosity are 
above the threshold L[z) becomes. 



n{y) = ^ [z{y)] . 



(35) 



dN ^ {d^f{l + z)dno'if{y)dy. 



(36) 




(37) 



while the total count yields. 



dN{z) 



(d^^dnoii + z)^ 



dy 
dz 



(38) 



dz 



The observed relativistic density per source in the universe will then be given by. 




(39) 



In general this function will also depend on 9 and ip. 



- 15 - 



which imphes that, 

P = Mg^, (40) 

if we assume the cosmological constant is zero and that all dark matter is associated with 
the sources. Otherwise, we would have to write pm — M-g^ , where pm represents only the 
matter associated with the sources themselves, not vacuum energy nor, for instance, matter 
distributed differently from the observed sources. These equations can be completely deter- 
mined in terms of the redshift once the function y = y(z) is found for a specific cosmological 
spacetime, allowing us to eliminate the affine parameter. 

For bolometric counts the lower limit L(z) is given by equation (30), while for counts 
in the bandwidth W we must use equation (31). If observations are in a limited bandwidth 
W, then the observed relativistic density per source of the universe must be written as, 

f^oi^) = EprwA^.w^rw(i+^)*r(^)^, (41) 

where 

Kiz) - / (42) 

^0 

and J^{z) is the completeness parameter determined by means of the bandwidth luminosity 
as given by equation (31), but with a luminosity function for each galaxy morphological 
type V and at different wavebands W . Notice that since the luminosity function conveys 
observational information about the number of sources at each redshift bin, it will effectively 
take into account the possible galaxy mergers taking place in those ranges, with the mor- 
phological population fraction P^{z) conveying the observational relative abundance of each 
morphological type in the same redshift ranges, and for each observing bandwidth W . In 
other words, number evolution and cosmological evolution are contained in the luminosity 
function, but number evolution will also be reflected in M.g{z). 

Remembering equation (23), the observed relativistic density per source can, for the best 
possible observations, be at most equal to its theoretical counterpart, iiq{z) < n{z). However, 
if we now suppose the existence of some form of dark matter, we should rewrite equation 
(23), since the total theoretical value of fj, should be given by the following expression, 

= ^ + (43) 

where fi^ accounts for the possible dark matter component. Therefore, from now on all 
theoretical quantities will be decomposed into an observational part, and a possibly unseen, 
differently distributed, component. 



-16- 



According to this definition, ji^ does not include dark matter galactic halos, as this is 
already included in the rest mass determination of each galaxy. So, /x^ may be baryonic 
intracluster dark matter, or non-baryonic dark matter, or vacuum energy, or still possibly, 
quintessence. 

In writing ^xq{z) we have also to consider that it is composed by summing up the 
individual contributions of the different bandwidths W , W, etc, where the observations are 
actually made. Therefore, we have that, 

^,{z)=f,^ + fif+lj.r + --- (44) 

This equation is valid provided that the different bandwidths W, W, etc, do not overlap, 
otherwise there will be double counting. However, in general we may expect significant 
overlap of galaxy counts over different bandwidths, meaning that equation (44) will convey 
incorrect results due to overcounting of the galaxies we observe. In order to avoid this 
situation, let us rewrite equation (44) as follows, 

Wf 

/^o(^) = Yl ^^i^)l^oi^)^ (45) 

w=l 

where 

a^{z) = 1, for W = l, (46) 

and 

a^{z) = b,„-{z) < 1, for W>1. (47) 

Here bn,{z) is the fraction of galaxies in waveband W > 1 which are not counted in wavebands 
1,2,...,(PF-1). 

For later usage it is convenient to define the luminosity density as being given by (Peebles 
1993, p. 120), 

j= / L<f>{e)de. (48) 

This is the average luminosity of the sources per unit of proper or comoving volume. 



2.3. Observed Average Densities 

The observed average density of the universe is defined as follows. 



(P) . (49) 
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where Vq is the total observed volume. As discussed in Ribeiro (2001, and references therein), 
equation (49) does not define an unique quantity, but at least three, as the observational 
volume can be constructed with any of the three distance definitions given by the reciprocity 
theorem (24). Consequently, we have that 

Vo = ^ndo^ (50) 

where may be any of the three distances, rf^, rf^, or d^.^ As at the observer's position 
C{q) the redshift is zero, these definitions must imply that 

= Po- (51) 

Notice that the observational quantities so far used, namely n, and (p, are not average 
densities, but local densities, defined at the source's rest frame. 

To determine (p) the key quantity to be found is A^. Bearing in mind the results 
previously derived, the number of observed sources above a certain luminosity L{z) is given 
by, 

No{z, >L)^Y1 = E /' + '^^o -^"^ "^y^ (52) 

w=l w=l 

while the overall number of sources of any luminosity is, 

Wf Wf 

/•y 

(53) 



N{z) = J2 awiz)N'''{z) = Y,aw I {dAf{l + z) dn^ dy. 

w=l w=l ''^ 



We can now write the observed and total average densities of the universe in terms of 
the previously defined quantities. The expressions yield, 

Wf y 

^P^o = 7^ E aw{z)J2Pv^i^)M.{z) r {d,)\l + z) d^o C dy, (54) 

and 

Wf y 

(P) = (P), + — ^ E «-(^) E r (d.)\l + z) dQo K dy, (55) 

47r(rfo) ^ Jo 



^ By means of the observed proportionality between redshift and distance, i.e., the Hubble law, in FLRW 
cosmologies one can also define the redshift distance, given by dz = cz/Hq, where c is the light speed and 
iJo is the Hubble constant. Therefore, a fourth distance definition is also available in the standard model, as 
well as another total observed volume. For non-FLRW metrics, one can only have a similarly defined dz if 
an equivalent to the Hubble parameter is also defined in the chosen geometry {e.g., Ribeiro 1994; Pompilio 
and Montuori 2002). 
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where (p)^ is the possible dark matter component. In redshift bins where varies httle with 
we have the following result, 

(p)o = J{P) - {P)a- (56) 

It is clear from equations (54) and (55) that building an average density from observa- 
tions is a very complex task which depends on many choices, like the appropriate distance 
definitions to be used, and a still to be determined function y — y{z) (if we do not assume a 
cosmology to begin with), which is certain to be non-linear. In addition, the complexity and 
high non-linearity of (p) shows the inherent difficulties of the task of observationally ascer- 
taining whether or not the universe is smooth when one follows a fully relativistic approach 
of this problem (Ribeiro 2001). 



2.4. Mass-to-Luminosity Ratio 

The mass-to-luminosity ratio at some redshift depth is defined as. 



M/L = ^ = Mgj, (57) 

while its observable counterpart is given by, 

[M/L], = M,t (58) 

If the average luminosity of the universe is similarly defined as the average density (p) , we 
have that, 

(i) = (59) 



where Lg is the average galaxy luminosity, given by 

V 

L,(^)=.^P,(^)L,(^). (60) 



v=l 



Here L„ is the typical luminosity value for galaxies of a certain morphological type v, and 
Py is the galactic fractional population class (see eqs. 13 and 14). 

From the equations above, the average mass-to-luminosity ratio at some redshift depth 
is given by. 
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3. Luminosity Function Data 

The discussion presented so far has focused on galaxy number counts, which, aside from 
the redshift, is one of the two observational quantities (observer area distance, or luminosity 
distance is the other) which determine the geometry and mass-energy distribution of the 
Universe (assuming spherical symmetry - if this is not assumed, then we need two other 
quantities, cosmological proper motions and null shear parameters - see Ellis et al. 1985). 
Thus, the ideal way to determine the mass-energy distribution and the geometry would be by 
solving the field equations directly with redshift, number count and observer area distance 
data. However, as already mentioned, our ignorance of luminosity and number evolution, 
and adequate model of which would have to be assumed, is an obstacle to implementing 
this in a compelling way. Because of this, and since the luminosity functions available in the 
literature presuppose a definite FLRW cosmology, the first and easiest sensible application of 
the general theory we have summarized above is to test luminosity functions for consistency 
with observational relationships determined by the cosmologies they have assumed. In other 
words, in any specific cosmological metric we need to make both sides of the equation 

[geometric cosmological model] ^ [parameterized data] 

consistent, as we will be dealing either with quantities which are determined only observa- 
tionally {e.g., 0), or with theoretical quantities derived from the chosen cosmological model 
{e.g., n, N), or both {e.g., dN/dz, M/L) (Stoeger 1987, p. 299). 

The luminosity function provides an excellent observational tool for testing cosmological 
models, as it conveys information on how the local number of galaxies n vary for different 
types of galaxies, its environmental dependence, how n changes with the cosmic epoch (see 
eq. 35), and also a possible indirect estimate of the amount of dark energy present. The 
luminosity functions also provide a statistical description of the galaxy population. The 
functional form of and its parameter values are extracted from magnitude-limited redshift 
surveys of galaxies by means of various statistical techniques (Peebles 1993, p. 119; Binney 
and Merrifield 1998, p. 162; Peacock 1999, §13.3). Here, however, we shall only be interested 
in collecting from the literature the shape, normalization, and the luminosity scale of for 
the purpose of using it in our equations for the mass-energy density and the relativistic mass- 
energy density per source, and checking that the results are consistent with the cosmology 
assumed in calculating the luminosity functions. Some of the information contained in the 
luminosity functions is tied to the galaxy morphological or spectral class (see below) , which 
means that at least a simplified morphological division must be assumed. 

As discussed by Chen et al. (2001), various deep redshift surveys have yielded consistent 
measurements of the luminosity function for galaxies at 2; ^ 0.75 (Marzke et al. 1998; Lin 
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et al. 1999; Cole et al. 2001), and this is precisely the redshift depth we are interested in 
here (see §1). Impey and Bothun (1997) have also compiled valuable observational data on 
0. Since these studies yield recent data on 0, they seem reliable for giving the luminosity 
functions for this redshift range. 

3.1. Morphological Types 

As it is well known, galaxies have different morphologies, which vary from ellipticals 
and spheroidals to spirals and irregularly shaped galaxies (Binney and Merrifield 1998, ch. 
4; Elmegreen 1998, ch. 2; Sparke and Gallagher 2000, §1.3). Redshift surveys dealing with 
thousands of galaxies tend, however, to automatically classify galaxies according to their 
spectral energy distributions (Lin et al. 1999, p. 539; Cole et al. 2001, §5.3), or by means 
of a simple three-type scheme, as E/SO, spiral, and irregular/peculiar galaxies (Marzke et 
al. 1998). It is important to mention that spectral type classifications are not necessarily 
straightforwardly related to the galaxy types of the Hubble diagram, and its extensions, and, 
since many surveys do adopt such a scheme, for the purposes of this paper it seems enough 
to follow Lin's et al. (1999) spectral morphological types. So, from now on we will assume 
V — ?> (see eq. 13), for three loose classification types: v — 1 stands for early type galaxies^ 
roughly E/SO; v = 2 accounts for intermediate type galaxies, something like Sa-Sb; and f = 3 
represents late type galaxies., which would be Sc-1 types. Notice, again, that this equivalence 
between the galactic spectral energy distributions and morphological types is not strict, and, 
therefore, should not be taken as giving a precise classification. 

The classification scheme above is not unique. For instance, with the data presented 
by Impey and Bothun (1997) one could think of only two galactic types: f = 1 for galaxies 
in clusters, which tend to be dominated by spirals, and u = 2 for field galaxies, which 
are usually ellipticals. One could even think of other classification schemes, hke luminous 
galaxies {v — 1), low surface brightness galaxies {v — 2), and luminous giants galaxies 
{v — 2)), or still V — 1 ior red galaxies, and blue galaxies as v — 2. Perhaps, combinations 
of all these schemes are possible, although in such a case the problem of overcounting could 
become more acute. 

Nevertheless, whatever classification scheme is chosen when applying the theory de- 
veloped above, the bottom line is that the adopted morphological types will be dependent 
basically on the sample, or on the sampling method, used in the chosen redshift survey. In 
most cases it will probably not involve strict identification with the galaxy morphologies of 
the Hubble diagram. So, from now on the term "morphological type" will have a broad 
meaning, referring to some sort of classification obtained by differentiating the galaxies con- 
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tained within a certain sample. Such a classification can be based either on spectra, or on 
the galactic shapes as given by the Hubble diagram, or on any other galactic observational 
feature considered important by some sampling method. 

3.2. Shape 

In finding the shape of the luminosity function, most, if not all, recent redshift surveys 
fit their data by using Schechter's (1976) elegant form, which summarized earlier findings by 
Zwicky and Abcll that the luminosity function is a broken power law, ~ L~^'^ at L < L^,, 
bending to ~ at the bright end (Peebles 1993, p. 119). Schechter's luminosity function 
is given by, 

(f){i) = 0, Te-^ (62) 

where 0* characterizes the space density of galaxies, and a is the asymptotic slope of the faint 
end of the luminosity function. The luminosity scale is given in £ (sec eq. 27), and tells us 
the luminosity above which galaxies are rare. These three constants are all observationally 
determined parameters. 

The power law index a for the faint end slope is part of the shape determination of 
the luminosity function, and, according to the latest observations, it does change with the 
morphological class of the surveyed galaxies, whether they are field galaxies, dwarfs, giants, 
spiral-rich galaxy clusters, luminous galaxies or low surface brightness galaxies (Impey and 
Bothun 1997). The analytical form of equation (62), however, does not appear to vary with 
galactic properties. Only the values of the parameters do. 

In the early days of the work on the luminosity function, all parameters were fitted to 
data only for local galaxies, where evolution is not considered to be important. However, as 
the surveys went deeper and deeper it became clear that possible luminosity evolution had 
to be taken into account. This can be done by allowing some of the fitted parameters to 
vary with the redshift (Lin et al. 1999), or by comparing their values at different redshift 
bins (Chen et al. 2001). Obviously, a great deal of work is needed to model luminosity and 
number evolution accurately, such as that pioneered by Totani, Yoshii, and their collaborators 
(Totani and Yoshii 2000; Totani et al. 2001). We arc not concerned in this paper with doing 
this. In order to illustrate how our framework can be used, we provisionally accept the 
rough evolutionary models assumed in the CN0C2 survey and test them for cosmological 
consistency (see below). 
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3.3. Normalization and Luminosity Scale 

Once the shape is determined, which includes determination of the faint-end slope, that 
is, the a parameter, we are still left with two unknowns, the normalization factor 0* and 
the luminosity scale L*. Various determinations of 0* show it not to be very sensitive to 
different FLRW models, even to those with non-zero cosmological constant, meaning that 
up to 2; ~ 1 these measurements appear to be robust. The same seems also true for a and 
L*. The major differences in these quantities seem to come from measurements of different 
morphological types and in different bandwidths (Binney and Merrifield 1998, p. 167; Lin et 
al. 1999; Cole et al. 2001). So, the present data seems to indicate only a weak dependence 
of these parameters on the cosmological models (Cross et al. 2001). 

Redshift surveys usually present the fitted L* parameter in terms of its characteristic 
absolute magnitude M^. So, in order to use these we have to rewrite equation (26) for the 
luminosity function as the number density of all galaxies between absolute magnitudes M 
and M + dM. For Schechter's form (62) that becomes (Binney and Merrifield 1998, p. 163), 

(t){e) d£ = 0, re-^d£ = (0.4 In 10) (f)^ioOA{i+a){M.-M) ^_-^^qOA{m.-m)^ ^ ^^^^ 

(63) 

since the definition of absolute magnitude implies that I — lO°-^(^* ^\ 

3.4. The CNOC2 Redshift Survey Data 

We have chosen to extract data to be used here from the Canadian Network for Ob- 
servational Cosmology Field Calaxy Redshift Survey (CN0C2), presently the largest such 
a sample at intermediate redshifts. The CN0C2 sample for luminosity function (LP) de- 
termination contains over 2000 galaxies in the range 0.12 < z < 0.55, and with apparent 
magnitude in the red band within 17.0 < Rc < 21.5. Their fitting parameters were confirmed 
at 2; ~ 0.75 (Lin et al. 1999; see also Chen et al. 2001). The authors also claim their fitting 
to be valid with little difference for < 2; < 1. CN0C2 data were obtained and fitted in 
three bandwidths, B^b, Rc, and U, and are available for FLRW models with qq = 0.5 and 
qo = 0.1. Therefore, by providing data in the redshift range we are interested in for differ- 
ent bandwidths, and presupposing the Einstein-de Sitter cosmological model, the CN0C2 
survey is the most suitable to exemplify how we can use this theoretical-observational frame- 
work for testing the consistency of the LP's with the assumed cosmology, for determining 
how much dark matter there is distributed differently than what is assumed associated with 
galaxies, and for eventually determining the range of validity of simplifying assumptions at 
low redshifts, among other things. 
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The required data for the CN0C2 LP's assuming EdS universe are summarized in the 
tables 1, 2, and 3, where v denotes the three morphological types in which the sample 
population is divided: early, intermediate and late spectral types (see §3.1 above), and 
the proportion of each galaxy population as compared to the whole sample. This is a fixed 
value in all CN0C2 redshift range. M'^ and Q are two constants of the linear equation for 
the characteristic luminosity scale evolution (Lin et al. 1999, §3.2), 

M,(^) =M:-Q(^-0.3), (64) 

and 0* is in units of h'^ Mpc"'^ mag~^. 

In computing these parameters of the LF from the CN0C2 survey, it is important to 
note that Lin et al. (1999) also adopted a model for galaxy number density evolution and 
assumed that the parameter a, which governed the shape of the LF at the faint end, does 
not change with redshift for a given morphological type. They were also careful to consider 
surface density selection effects, and estimate the incompleteness of their sample, especially 
at the higher end of their redshift range. Finally, they performed various consistency checks 
to make sure that the LF parameters derived adequately represent the data for the redshifts 
indicated. It is reassuring that Cohen (2002) finds that her results for the Caltech Faint 
Galaxy Redshift Survey agree very well with the CN0C2 results, including for the faint-end 
slope of the LF's. Similarly, the results from the Deep Groth Strip Survey (Im et al. 2002) 
are in agreement with the CN0C2 results, especially for the luminosity evolution parameter 
Q. The one discrepancy that has emerged is that the significant number evolution for late- 
type galaxies in the CN0C2 results is not found by either Cohen (2002) or Im et al. (2002). 
Im et al. (2002) speculate that this may be due to the fact that the CN0C2 samples were 
selected purely by colors. This means at low redshifts {z < 0.4), where the difference in 
colors for various galaxy types is small, photometric errors can easily bump blue galaxies 
into the red-galaxy samples. Other than this, the CN0C2 LF results seem well supported 
for this redshift range. 

The determination of the LF's with CN0C2 data employed the comoving volume instead 
of the proper volume (Huan Lin 2002, private communication). Therefore, to be usable for 
us here we must apply an appropriate volume conversion factor to the results listed in tables 
1-3 (see §4.3 below). 

4. CNOC2 and Einstein-de Sitter 

As already mentioned, we now want to use our observational-theoretical framework first 
to test the consistency of the CN0C2 LF results with the cosmological model it assumes. 
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which is EdS. This has the added advantage that EdS is mathematically the simplest of 
the FLRW models. Furthermore, it is a frequently used case in the LF's available in the 
literature. The extension of this type of analysis to other more general FLRW cosmologies, 
including those with a non-zero cosmological constant, are the subject of a forthcoming paper 
(Stoeger and Ribeiro 2003). 



4.1. The Metric 



Let us write the EdS metric as follows (c — G—l), 

ds'' = dt'' - a\t) [dr'' + {dO'' + sin^ Odip'')] , 
where a{t) is the scale factor, given by 

Hq is the Hubble constant, and the local density is 

1 

The solution of the past light cone equation, dt/dy = —a{t)dr/dy, may be written as. 



t{y) + 



5Ho 



1/3 



1/3 



(65) 



(66) 



(67) 



(68) 



which allows us to write equation (66) along the null cone parameterized by r, as follows, 

a[t{r)] = - i^—j -r . (69) 

Detailed calculations on this spacetime along the lines followed in this article are widely 
available (see, for instance, Ribeiro 1995, 2001). Therefore, here we will only summarize 
some important results required in what follows. 



4.2. Observables 



Equations (66) and (69) allow us to obtain the function z = z[r{y)], that is, an expression 
for the redshift where the affine parameter becomes implicit. In EdS it may be written as. 



1 + z 



ait) 



18 
Ho) 



18 



1/3 



(70) 
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Applying equation (25) to the metric (65) gives us the area distance in this cosmology, which, 
by means of equation (70), yields. 



— ra{t) 



VTT^ - 1 
(1 + ^) 



3/2 



(71) 



since along the past null cone the scale factor may be written as. 



Hn 



2/3 



9{l + z) 



In EdS geometry, for comoving sources, — 6q, we have that, 

-k'^Qoh = —k^ = —dt/dy = a{dr/dy), 



(72) 



(73) 



which substituted into equation (1), and also remembering equation (16), gives us the total 
relativistic density per source (17) in the Einstein-de Sitter cosmology, as written below. 



dN 3Hn 



(d^) dQo dz Stt 



(74) 



From these expressions, the results below readily follow. 



dN 
dz 

N(z) 



3dn, (/TT^-i)' 



dQn 



ttHqMc 



1 



Stt 



(75) 



(76) 



(77) 



Finally, with either equations (18) and (74). or (11) and (75), wc can calculate the relation- 
ship between the affine parameter and the rcdshift in this cosmology. It may be written as 
follows, 

(78) 



dy^ 1__ 

dz Ho(l + z) 



7/2 ■ 



4.3. Volumes 



Before we proceed with the tests, there is one more issue to be taken care of: volume 
definitions. This is a tricky issue in relativistic cosmology as in general relativity volumes are 
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not uniquely definable quantities, and they may vary quite substantially from one another 
even at ^ 0.1 (Ribeiro 2001). Therefore, to avoid confusion, let us state clearly what 
volume definitions we have been using. 

Equation (35) states that both the luminosity function and number density are defined 
using local volumes, that is, volume elements defined at the rest-frame of a galaxy at a given 
redshift z. Prom equations (11) and (12) it is clear that this volume is given by the following 
expression, 

(^Mocai = dl.dSo = (1 + z)dy{d^fdno- (79) 
Considering equation (78) it follows that in EdS cosmology this reduces to, 

dz 

C^^ocal = {d^fdno— ^- (80) 

Ho{l + z) 

Sometimes observers refer to their data being reduced against the real, or proper volume, 
obtained by means of the spatial portion of metric (65), as follows, 

dVpj-oper — CL^r^dr sin 9 d9d(p. (81) 

Remembering equations (70) and (71), one can easily show that proper and local volumes, 
that is, equations (80) and (81), are equal. Nevertheless, we believe that since both the 
luminosity function and number density are defined locally, it makes more conceptual sense 
to refer to the local volume, as given by equations (79) and (80). In addition, it is important 
to stress that the local and proper volumes are in general not equal to the observed volumes 
defined by equation (50), as observed volumes use the three observational distances given by 
the reciprocity theorem (24). 

When carrying out their data reduction, observers also used to refer to a survey volume, 
defined as (Sandage et al. 1979), 

14urvey=^^if^o(^y, (82) 

where Wmax is the maximum velocity obtained in the survey, being calculated by means of the 
Doppler velocity-rcdshift approximation U = cz, valid only for small z. As it is well known, 
it is only under this approximation that we can write the rcdshift-distancc law, that is, the 
Hubble law, as a velocity-distance law (Harrison 1993), since the velocity-distance law is a 
general result of an isotropic homogeneous expanding universe, whereas the redshift-distance 
law is a mere approximation of the form, valid at small redshift. That, in turn allows us to 
obtain the survey volume. Such a procedure not only employs the low redshift approximation, 
but also implicitly assumes the redshift distance d^ as the distance definition used in V^urvey 
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(see footnote at page 17). Therefore, if one follows this path, model independent observable 
distances, like or dj^, are related to d^ only in an approximate manner in EdS cosmological 
model, dz approximately scales only with d^, in addition to following the same asymptotic 
behavior (Ribeiro 2001). 

Finally and most importantly, nowadays observers mostly use the comoving volume 
when determining LF parameters of galaxy redshift survey data. In this paper we have been 
using the proper volume in all calculations because it is more conceptually natural. Thus, we 
require a conversion factor to fit the CN0C2 data into our equations. Defining a comoving 
volume requires a metric function in comoving coordinates. In EdS cosmology it is obvious 
from equations (72) and (81) that, 

4 _3 

^^^roper — „ 2 ("^ ^^omoving) (^^) 

which shows clearly the conversion factor necessary in our calculations. 



4.4. Consistency Equations 

Before we write down the consistency equations for the LF parameters with respect to 
the underlying cosmology they assume, let us first bricfiy summarize the general formalism 
adopted so far. Let T be some theoretical quantity completely determined by theory, and 
To its measurement. Then, if J is the completeness parameter, solely written in terms of the 
LF, we have that 

To = JT. (84) 
The observational quantity is obtained at each bandwidth W. So, 

To = J2(^wT^: (85) 

w 

where is a quantity directly built from observations for instance). 

The dark matter is assumed to be in two components: (i) galactic dark matter, ac- 
counted for by J, as it ought to be included when measuring A4g, and (ii) intergalactic, 
possibly non-baryonic, dark matter, accounted for by the term T^. So, does not include 
any kind of galactic dark matter. In addition, if there is intergalactic dark matter, it must 
be added to Tq in order make up for the whole energy-matter as given by T. Consequently, 
equation (84) must be rewritten, 

T = T, + J2a^j^. (86) 

w 

Bearing those points in mind, we can now proceed in building our consistency equations. 
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4.5. First Consistency Equation: Local Density and Number Evolution 

4.5.1. Proper Volume 

Equations (16), (35), and (77) will give the theoretical behavior of the luminosity func- 
tion in EdS cosmology. In proper volume units, and in correct dimensions, it may be written 
as, 

*^» = i^(^ + ^'° («^' 

Remembering equations (13), (42), (46), and (47), and that the luminosity function is deter- 
mined for each galaxy morphological, or spectral, type, and at each observational bandwidth, 
the observational counterpart of the equation above is given by, 

i^M,), = J2a^{r'M,) = X]a^^P-(;.)A<,(^)C(^)- (88) 

w=l w=l v=l 

These two equations provide the first consistency test between Einstein-de Sitter theoretical 
prediction (eq. 87) and observations (eq. 88) given by the LF. They basically probe the 
behavior of the local density and number evolution of sources at higher redshifts. Com- 
paring them also gives us an estimate of the intergalactic dark matter component. In fact, 
remembering equation (86), its dependence on the redshift can be rewritten, 

2 Wf V 

Pa = + - E E p:riz)M.iz)Kiz)- (89) 

W=l v=l 

Notice that for other more general FLRW cosmologies, equation (87) will have additional 
terms due to the deceleration parameter qq and a possibly non-zero cosmological constant, 
while equation (88) will remain the same, inasmuch as it is built from observations. 



4.5.2. Comoving Volume 

For a comoving volume, we require the conversion factor (83) to be applied to equation 
(87). In correct dimensions, it yields. 

This is the quantity to be compared with the observational results as given by the calculation 
of equation (88) from the CN0C2 data. 
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4.6. Second Consistency Equation: Differential Number Counting 

Considering equations (17), (38), (78), and remembering equation (85), we can write the 
observational relativistic density per source, which is proportional to the differential number 
counts, of an Einstein-de Sitter universe in the volume units in which the LF is written, as 
follows, 

^^o{z) = — -^-^E E Pv^{z)M.{z)^:{z). (91) 

-"oil + w=l v=l 

This expression together with equation (74) provides the second consistency test in EdS 
cosmology. 

It is important to mention that equation (91) evaluates fiQ by means of the LF data, 
whereas if data were to be provided on differential number counts, /xq could be calculated 
by means of equation (22), where fewer hypotheses are "plugged" in for the observational 
data. This is a point we want to stress. Equation (22) only gives differential number counts 
for either some specific angular regions or over the whole sky, summed over all wavebands, 
which means that there is no need to evaluate the proper or comoving volume, as it is the 
case for the LF. Therefore, the best prospect for observational cosmology tests which use 
redshift survey data lies in the determination of the relativistic density per source of the 
universe, which in turn requires differential number counting observations as free as possible 
from any presupposed cosmological model. 

As in §4.5 above, the dark matter component associated with the relativistic density 
per source of the universe may, in correct dimensions, be written as, 

= - TTTT^^E E P.^{^)M.{^)K{^)- (92) 

Ho[l + z) 

It is also worth mentioning that equations (54) and (58) provide ways of measuring the 
average density and the mass-to-luminosity ratio, respectively. However, from a theoretical 
viewpoint they are just extensions of the two tests above, with no new quantities involved. 
Therefore, we shall leave their testing to future work. 

4-6.1. Comoving Volume 

As with the local density, to obtain the theoretical relativistic density per source of the 
Universe in comoving volume units we need to apply the conversion factor (83) to its proper 
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volume counterpart. Prom equations (74) and (83), it is clear that the theoretical prediction 
of the model for /i (comoving) may be written as, 



3cHq 
SttG 



vTT 



9(1 + ^)' 



cHq 



1 + z) 



-5/2 



(93) 



whereas the observational quantity is just equation (91), since Lin et al. (1999) CN0C2 LF 
parameters are already written in comoving volume units. 



5. Consistency Tests with the CNOC2 Data 

Before we start to actually test the CN0C2 data against the predictions of the EdS 
model, we still need to describe the means of evaluating equations (28), (29), and (33) in 
some waveband W . The critical equation in this respect is the waveband limited version of 
the selection function (eq. 28), since the other two come from this one. Its waveband version 
may be written as follows, 

/•oo 

^"^{z) = / (t>'^{^)dL (94) 

J Lw{z)/L,, 

Here L-w{z) is the minimum absolute luminosity in the filter waveband W , as determined 
at the observer, of sources detected in the given sample. The source spectral function, 
J [z/(l + z)], the proportion of radiation it emits in each frequency at the source itself, is 
also observed here, but corrected by the factor {1 + z) to account for its cosmological redshift. 
By examining equation (31) one can only conclude that Lw{z) must be written as follows, 

/■oo 

L^{z)= L{z)W{u)J[u{l + z)]du. (95) 

The luminosity evolution of sources is represented by their intrinsic luminosity as a function 
of redshift, L{z), as given in its rest frame, multiplied by J^[j^{l + z)]. But, without some 
theory for intrinsic luminosity evolution, the observer can only have access to the product 
{L{z)J' [1/(1 + z)]} as seen through the filter W. In principle, J" is also a function of the 
morphology and age of the source, and this is reflected in samples which make distinctions 
on the basis of source morphology and spectral type. 

Equation (31) may be rewritten in terms of absolute magnitudes Mw{z) if we consider 
the expression above, and the reciprocity theorem (24), yielding. 



(96) 
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where Cw is the zero-point-magnitude-scale constant relative to the filter W. This is a general 
equation, valid for any cosmological model. The (1 + z) factor appearing in the denominator 
accounts for the rclativistic correction, as compared to the Euclidean case. Schechter's LF 
(62) together with equation (96) allows us to numerically compute the selection function by 
means of simple quadratures of equation (94). 

Alternatively, the selection function can also be evaluated if we consider equation (63) 
for Schechter's LF written in terms of the absolute magnitude. Thus, we can rewrite the 
VT-filter selection function, as follows, 

ilj'^{z)= (O.41nlO)0aO°-^(^+")(A'^*-Mexp -10°-<^*-^^) dM^. (97) 

J —oo 

Considering equation (64), which models the luminosity evolution in the CN0C2 survey, 
this equation takes the form. 



ip'^iz) = 0.4 In 10 0, (98) 

"Mwiz) 

X ' 

-oo 



-^Q0A{l+a)[Mi-Q{z-0.3)-Mw] (-^Q0A[Mi-Q{z-0.3)-Mw] \ dJf^ 

-oo ^ 



In either case 0*, M^, Q, a will take the fitted values as given in tables 1-3, the absolute 
magnitude range for each filter used in the LF derivations of the CN0C2 sample, for Hubble 
constant h ^ 1, are, -23.0 < < -17.0, -22.0 < < -16.0, -22.0 < M^^ < -16.0 
(Lin et al. 1999, tables 1-2), and the magnitude zero-point constant of equation (96) 
in each filter are, C^^ — —13.64, Cu — —14.06, Cab — —48.60 (Fukugita, Shimasaku and 
Ichikawa 1995). 

Once t/j is known by either method above (we used the second method, with absolute 
magnitudes, as it requires fewer data manipulations), we can calculate the summation term 
for all filters and all galaxy morphologies, 

D(z)^ h^Y.P^M.i^'^iz). (99) 

w=l v=l 

(Rc,U,AB) 

Here fe^^ = = b^B = 1/3 (see eq. 47), since the same CN0C2 sample was used at 
different bandwidth observational windows for obtaining the LF parameters. The dynamical 
masses for each morphological type are Mi = 0.5 x IO^^Mq, M.2 = 0.3 x IO^^Mq, M3 = 
0.1 x IO^^Mq (Sparke and Gallagher 2000, p. 204, 264; see also §2.1.2 above). The results 
are summarized in tables 4-7. 
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5.1. First Test 

The first consistency test compares the theoretical predictions of the local density in 
comoving coordinates (eq. 90), and its observational equivalent. They are respectively given 
by, 

mMg = — ^ = 37.8 X lO^A^o Mpc-2/i^ (100) 
dttG 

and (see eq. 88), 

{^M,\ = D. (101) 

Results of the D-term are found in table 7. One can clearly see that the observational results 

arc approximately constant around 0.5, which gives the luminous matter up to 2; = 1, as 
detected, and extrapolated, by the CN0C2 galaxy redshift survey, as being about 1.3% of 
the critical mass of EdS universe model. Since this quantity is approximately constant, 
one can conclude that the LF determination in this catalogue is, at least in this respect, 
consistent with the theoretical predictions. In addition, this percentage of luminous matter 
as compared to the total mass seems to fall in line with some of the current expectations 
that the Universe is filled with large quantities of dark mass, much of it non-baryonic, and 
not distributed as the galaxies are. 

One can also see a gradual rise in D which might have some marginal significance. As 
equation (100) indicates, theoretically ^M.g should be constant. But, observationally, from 
the LF parameters it is not because Lin et al. (1999) have included some number density 
evolution, which is reflected in (that is clear from the discussion in their paper). Thus, 
is a function of redshift, and increases slightly with redshift. However, since we arc holding 
M.g constant, their product increases with redshift. In order to keep the product constant 
we would have to give Jvig a dependence on redshift in order to cancel the dependence of ^0 
on it. 



5.2. Second Test 

The second test involves comparing the theoretical expression for the relativistic density 

per source in a comoving volume, which is equivalent to number counts, as given by equation 
(93), and its observationally-derived-LF expression (eq. 91), which may also be written as. 

The functional relationship of the observational and theoretical results against the redshift 
are respectively given in the third and fourth rows of table 7, where one can clearly see a 
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similar downward trend for both quantities at higher redshifts. Nevertheless, the theoretical 
predictions do tend to decrease more quickly than the observational results. At the range 
0.1 ^ z ^ 0.25 the observational results tend to level off, restarting a decrease only at 0.3 
(sec figure 2). Since /io gives information about the differential number counts (see eqs. 17 
and 22), it seems thus reasonable to conclude that the theory does not adequately account 
for the galaxy counts at higher redshifts. This may reflect either the inadequacy of the EdS 
model or of the number evolution model Lin et al. (1999) are employing in the calculation 
of their LP's. 



6. Conclusions 

In this paper wc have discussed the theory connecting the rclativistic cosmology number 
count theory with the astronomical data, practice and theory behind the galaxy luminosity 
function (LF). We started from EUis' (1971) general, and model independent, relativistic 
expression for number counts at a point along the observer's null cone and derived expressions 
for the relativistic density per source of the universe, differential number counts, mass- 
to-luminosity ratio, and local and average densities. We then linked those expressions to 
the current definition adopted in observational cosmology for the luminosity and selection 
functions. The equations were constructed in such a way as to reflect the fact that the 
luminosity function is determined only within certain bandwidth ranges, and for some galaxy 
morphological types. 

We then discussed the current determination of the LF parameters in order to connect 
the theory to current practice. The resulting theoretical-observational framework for con- 
fronting theory with cosmologically relevant observational results can be used in a number of 
ways. We can factor the cosmology out of the LF's to obtain the differential number counts; 
we can investigate the range of validity of various low-redshift and Euclidean approxima- 
tions; as we have done here, we can take LF results and quickly check their consistency 
with the cosmology they assume for a given range of redshift; we can use LF number counts 
results, with or without assumed models of luminosity and number density evolution, and 
either fit them to a given cosmology or calculate the cosmologies they determine. In order 
to accomplish the latter we also need observer area distance (or luminosity distance) vs. red- 
shift data. As we mentioned in passing, an analogous theoretical-observational framework 
can be easily constructed for this key cosmological parameter. In this paper, as a simple 
example of how the general theoretical-observational framework can be used, we examined 
the consistency between the LF parameters of a redshift survey and some of the key equa- 
tions of the cosmology they assume. We chose the CN0C2 intermediate redshift survey, and 
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the EdS determined parameters that Lin et al. (1999) calculated from it, to illustrate this in 
a simple way. After specializing our general model to EdS, we carried out two consistency 
tests, one related to the local density and number evolution, and another dealing with the 
differential number counting, comparing the theoretical predictions with the observational 
results derived from the LP determination of the CN0C2 survey. The results show a gen- 
eral agreement between theory and observation in both tests. Nevertheless, the second test 
indicated a possible excess of galaxies at 0.1 < 2; < 0.4, which either means that the EdS 
cosmology does not accurately represent the data, or that the evolution model assumed by 
Lin et al. (1999) is not adequate. 

It is important to point out that the actual observations used in the CN0C2 survey 
LP determination were carried out only within the range 0.12 < z < 0.55, which means 
that in most of its observational redshift interval there appears to be an excess of galaxies 
as compared to EdS predictions. Por z > 0.4 our results indicate that this excess tends to 
disappear. In order to account for this extra mass we can either assume that the universe 
has a closed PLRW geometry, or hypothesize, as many others have been doing recently, 
a non-zero cosmological constant as a possible extra source of mass-energy, or modify the 
model of number evolution Lin et al. (1999) are employing. 

The fourth option, a possible inadequacy of the EdS cosmological model, appears to us 
less likely, since the LP determinations seem to have only a weak dependence on the chosen 
model. If we were to be suspicious about the inadequacy of the cosmological model, due to 
this weak dependency we should perhaps call into question not the EdS universe, but the 
whole set of standard PLRW cosmologies. It seems to us that our results do not support 
such a more radical viewpoint, but this is something which one always need to bear in mind. 

We are grateful to Chris Impey for illuminating discussions on the methodologies behind 
data reduction and analysis of redshift surveys, to Huan Lin for clarifications regarding 
the CN0C2 survey, and a referee for useful remarks. MBR wishes to thank the Vatican 
Observatory Research Group for their kind hospitality while performing this research, and 
to acknowledge the financial support from Brazil's CAPES Poundation. 
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] C(q): point q of observation 
along C 



i C (observer' s world line) 

Fig. 1. — Section of a bundle of light rays subtending a solid angle dQo as seen by the observer 
in C{q). The affine parameter displacement dy corresponds to a local distance variation of 
{—k"'Ua) dy in the rest-frame of the galaxy at a point P down the null cone (EUis 1971). 
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Fig. 2. — This plot shows the theoretical and observational results for the relativistic density 
of the universe (y-axis) vs. the redshift. The lower curve, made of points joined by dots, 
shows the observational results of //q, derived from the CN0C2 redshift survey, whereas the 
full upper curve is the theoretical prediction for fi{z) in EdS cosmology. One can clearly 
see that although fiQ has the same downward tendency as /i, the theoretical results decrease 
more rapidly in the range 0.1 ^ 2; ^ 0.4. For z > 0.4 an equivalent-to-the-theoretical- values 
downward trend is resumed. At z = 0.05 /iq is about 1.3% of fi, becoming 1.6% at z = 0.25, 
only to return to be 1.3% oi n at z = 0.4 — 0.55. Notice that the CN0C2 sample used 
in the LF parameters determination went only up to 2; = 0.55, and the LF results were 
considered valid in the range < 2; < 1 by means of extrapolations (see §3.4 above, and 
Lin et al. 1999). These results indicate a broad qualitative agreement between theory and 
observations, apart from the region 0.1^2:^0.4, which may imply that either the chosen EdS 
cosmology provides an inadequate representation of the observational differential number 
counting data, or the CN0C2 survey detected an excess number of galaxies in most of its 
nominal observed redshift bin of 0.12 < 2; < 0.55 as compared to EdS predictions. 
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Table 1. LF fit in CN0C2 sample B^^ band (go = 0.5). 



V a (f)^ Q Ml P, 



1 0.08 0.0203 1.58 -19.06 0.29 

2 -0.53 0.0090 0.90 -19.38 0.24 

3 -1.23 0.0072 0.18 -19.26 0.47 



-41 - 



Table 2. LF fit in CN0C2 sample Rc band (go = 0.5). 



V a (f)^ Q Ml P, 



1 -0.07 0.0185 1.24 -20.50 0.29 

2 -0.61 0.0080 0.69 -20.47 0.24 

3 -1.34 0.0056 0.11 -20.11 0.47 
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Table 3. LF fit in CN0C2 sample U band (go = 0.5). 



V a (f)^ Q Ml P, 



1 0.14 0.0213 1.85 -18.54 0.29 

2 -0.51 0.0092 0.97 -19.27 0.24 

3 -1.14 0.0095 0.51 -19.32 0.47 
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Table 4. CN0C2 selection function^ results in the Rc band vs. redshift.^ 



z 








0.05 


0.0181 


0.0120 


0.0220 


0.12 


0.0182 


0.0121 


0.0221 


0.25 


0.0183 


0.0123 


0.0223 


0.4 


0.0185 


0.0125 


0.0225 


0.55 


0.0186 


0.0127 


0.0227 


0.75 


0.0187 


0.0130 


0.0230 


0.9 


0.0187 


0.0132 


0.0232 


1.0 


0.0187 


0.0133 


0.0234 



^Subscript numbers denote 
the spectral type morphology 
adopted for the CN0C2 survey 
(see §3.1). Units are Mpc~^. 

^^Redshift values were taken to 
agree with the nominal redshift 
interval and bins of the survey, 
and the claimed validity of the 
fitted LF parameters (see §3.4 
above, and Lin et al. 1999, p. 
550). 
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Table 5. CN0C2 selection function results in the U band vs. redshift. 



z 








0.05 


0.0180 


0.0119 


0.0284 


0.12 


0.0182 


0.0121 


0.0289 


0.25 


0.0186 


0.0123 


0.0297 


0.4 


0.0189 


0.0126 


0.0307 


0.55 


0.0192 


0.0128 


0.0317 


0.75 


0.0194 


0.0131 


0.0330 


0.9 


0.0195 


0.0134 


0.0341 


1.0 


0.0196 


0.0135 


0.0348 
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Table 6. CN0C2 selection function results in the Bab band vs. redshift. 



z 








0.05 


0.0182 


0.0122 


0.0252 


0.12 


0.0183 


0.0123 


0.0254 


0.25 


0.0185 


0.0125 


0.0257 


0.4 


0.0187 


0.0128 


0.0260 


0.55 


0.0189 


0.0130 


0.0264 


0.75 


0.0190 


0.0133 


0.0268 


0.9 


0.0191 


0.0135 


0.0272 


1.0 


0.0191 


0.0136 


0.0274 
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Table 7. D-term* sum (eq. 99) results vs. redshift. 



z 


D 


(i + ^)-'/'d 


(1 + ^)-'/''^ 


0.05 


0.47 


0.42 


0.89 


0.12 


0.47 


0.35 


0.75 


0.25 


0.48 


0.35 


0.72 


0.4 


0.49 


0.21 


0.43 


0.55 


0.49 


0.16 


0.33 


0.75 


0.50 


0.13 


0.25 


0.9 


0.51 


0.10 


0.20 


1.0 


0.51 


0.09 


0.18 



^Units are 10^ A^© Mpc'^. 

^Included for comparison (see eq. 93). 



